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Description 

FIELD OF THE INVENTION 

5 [0001] This invention relates to a device for classifying heart rhythms of a patient based on the shape (morphology) 
of the endocardial or epicardial electrogram. More particularly, a number of features (metrics) of the electrogram are 
measured to form a description of the shape of individual electrograms, and the metrics are then converted into a 
cardiac rhythm diagnosis by means of a pattern classification technique. 

10 BACKGROUND OF THE INVENTION 

[0002] There are a number of ways to discriminate different types of heart rhythm in an implantable cardiverter 
defibrillator (ICD). The most obvious way to do this is to use heart rate information alone, in the manner shown for 
example in U.S. Patent No. 4,475,551 to Langer et al. t entitled "Arrhythmia Detection and Defibrillation System and 
*5 Method", which issued October 9, 1984. The purpose of discriminating different rhythms in an ICD device is to allow 
the device to respond appropriately to each of them. Problems occur if there are two or more rhythms which require a 
different response from the device but cannot easily be discriminated from one another. 

[0003] The rhythms which cause most problems in existing ICD devices are atrial fibrillation (AF) and sinus tachy- 
cardia (ST). Both of these rhythms can result in a ventricular rate that is high enough for it to be considered to be 
20 ventricular tachycardia (VT). Thus, VT cannot be discriminated from AF or ST on the basis of ventricular rate alone. 
Additional information must be used to correctly identify these rhythms. 

[0004] As discussed in an article by Arzbaecher et al., entitled "Automatic Tachycardia Recognition", appearing in 
PACE , May-June 1984, Volume 7 (II), pages 541-547, and in an article by Jenkins et al., entitled "Tachycardia Detection 
in Implantable Antitachycardia Devices", appearing in PACE , Volume 7 (II), pages 1273-1277, November-December 
25 1 984, atrial fibrillation can be readily identified by looking at the ratio of the atrial rate to the ventricular rate. This requires 
the use of an atrial lead, however, which may be undesirable or unavailable. Also, this method is unsuited to discrim- 
ination of ST and VT, since both rhythms may exhibit a 1:1 ratio of atrial to ventricular rate. 

[0005] Timing methods have been proposed for distinguishing ST or AF from VT (see, e.g., U.S. Patent No. 4,830,006 
to Haluska et al. t entitled "Implantable Cardiac Stimulator for Detection and Treatment of Ventricular Arrhythmias", 
30 which issued on May 1 6, 1 989, and see, e.g., an article by Camm et al., entitled "Tachycardia Recognition by Implantable 
Electronic Devices", appearing in PACE , September-October 1987, volume 10, pages 1175-1190). These methods are 
of dubious accuracy and some of them also require the use of an atrial lead. 

[0006] The shape or morphology of the intracardiac electrogram can be used to achieve the desired discrimination 
between ST and VT. 

35 [0007] Template methods of morphology analysis provide a scalar measure of the difference between two shapes. 
An example of normal morphology is established, and this is compared to each heartbeat for analysis. A scalar bound 
or threshold on the normal variation in this scalar measure is established, and any measure that exceeds this threshold 
is considered to represent an abnormal beat. 

[0008] A number of template methods have been heretofore described (see, e.g., U.S. Patent No. 5,000,189 to 
40 Throne et al., entitled "Method and System for Monitoring Electrocardiographic Signals and Detecting a Pathological 
Cardiac Arrhythmia Such as Ventricular Tachycardia", which issued on March 19, 1991, and see, e.g., an article by 
Throne et al, entitled "A Comparison of Four New Time-Domain Techniques for Discriminating Monomorphic Ventricu- 
lar Tachycardia from Sinus Rhythm Using Ventricular Waveform Morphology", appearing in IEEE Tr. on BioMed. Eng ., 
Vol. 38 (6), at pages 561-570, (June 1991)). Template methods in general require high sampling rates and high process- 
es ing overheads. For these reasons they are unsuited for use in an ICD. 

[0009] Many attempts have been made to find a single metric that can be used for rhythm discrimination purposes 
(see, e.g., an article by Pannizzo et al., entitled "Automatic Methods for Detection of Tachyarrhythmias by Antitachy- 
cardia Devices", appearing in the Journal of American College of Cardiology at Vol. 11 , pages 308-316, February 1988). 
A probability density function (PDF) algorithm has been of some utility in diagnosing ventricular fibrillation (VF), but 
so has not proven successful for discriminating ST from VT. The PDF algorithm is discussed in the aforementioned U.S. 
Patent No. 4,475,551 to Langer et al., and in an article by Bach et al., entitled "Engineering Aspects of Implantable 
Defibrillators," appearing in a book edited by Saksena et al., entitled "Electrical Therapy for Cardiac Arrhythmias: Pac- 
ing, Antitachycardia Devices, Catheter Ablation", published by W. B. Saunders, Philadelphia, 1990, at pages 375-376. 
[0010] U.S. Patent No. 5,086,772 to Lamard et al. describes a method for combining two simple morphological 
55 features with timing information to improve the rhythm discrimination process. This method considers first the rate, and 
uses the morphological information only in a specified rate band, to classify individual heartbeats. The cardiac rhythm 
is then diagnosed on the basis of the classification of a number of successive heartbeats. The present invention is a 
generalization of the Larnard et al. morphological method but its application is not limited thereby. 
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[001 1] EP 0506230 describes a device for monitoring waveform complexes in which a comparison is made between 
a measured and a stored complex and a scoring method is used to assess the similarity of the morphology. 
[0012] It is, therefore, a primary object of this invention to provide an improved device for classifying heart rhythms 
by means of electrogram morphology. 

5 [0013] It is another object of this invention to provide an improved device for classifying heart rhythms by means of 
electrogram morphology, which method is usable in an implantable cardioverter defibrillator. 
[0014] It is a still further object of this invention to provide a device for classifying heart rhythms by means of elec- 
trogram morphology, which device employs a low sampling rate, requires very few multiplications per heartbeat clas- 
sification, consumes very little power in the classification process, and can be implemented both in digital and analog 

10 electronics. 

[0015] Further objects, features and advantages of the invention will become apparent as the following description 
proceeds. 



SUMMARY OF THE INVENTION 

15 

[0016] According to the present invention there is provided a device for classifying heart rhythms of a patient in 
accordance with the claims. 

[0017] Briefly stated, a number of metrics are used simultaneously to improve the accuracy of the classification 
decision process. Bounds on the normal variation of the metrics are established. These bounds define a decision 

20 surface which divides the measurement space into two regions. Measurement vectors which fall in one of these regions 
are considered to represent normal heartbeats, while those that fall in the other region are considered to represent 
abnormal heartbeats. Methods for achieving this have been discussed by K. Fukunaga in a book entitled "Introduction 
to Statistical Pattern Recognition", published by Academic Press, Inc., London, Second Edition, pp. 124-169, 1990, 
and by Duda et al., in a book entitled "Pattern Classification and Scene Analysis", published by John Wiley & Sons, 

25 New York, pp. 24-31, 1973. The particular metrics disclosed in the present invention require only simple operations 
and low sample rates to calculate. Further, the method for reaching a rhythm classification based on these metrics 
imposes a very small computational burden. Thus, this method is ideally suited for use in ICDs. 
[0018] The shape of the electrogram is quantified by measuring particular features thereof, such as amplitude or 
width. The value of these features are the metrics, and a set of metrics form a quantitative description of the shape of 

30 an electrogram. The shape of an electrogram is then indirectly compared to that of a known normal electrogram by 
comparing the metrics, and a given heartbeat is then classified as either normal or abnormal. VT may then be diagnosed 
if a number of successive heartbeats are abnormal. 

[0019] The metrics that are discussed here are simple ones that can be implemented either in analog electronics, 
or with low computional overhead on a microprocessor. They require only such techniques as filtering, peak-picking, 

35 thresholding, and integrating. Further, by use of a peak-holding sampling process, they can be calculated from a very 
low data rate digital signal. The rate of 80 Hz has been shown to be sufficient. A suitable set of metrics that may be 
employed in connection with diagnosing a rhythm are the area under the electrogram (hereinafter defined more com- 
pletely and referred to as "MOD"), the peak-to-peak amplitude (hereinafter defined more completely and referred to 
as "PKP"), and the amplitude of the largest negative peak (hereinafter defined more completely and referred to as 

40 "NVE"). 

[0020] The electrogram signal is filtered, preferably with a passband of two to thirty Hz (- 3dB). The lowpass section 
of the filter serves as an antialiasing filter, and smooths the transitions in the electrogram so that peaks may easily be 
identified. The highpass section of the filter serves to remove any baseline wander from the signal, so that a separate 
DC removal stage is not required. 

« [0021] The filtered signal is sampled at a rate that can be as low as 80 Hz, using a known peak-holding method, 
such as the one discussed at pages 217-220 in the book The art of Electronics" (2nd Ed.), authored by P. Horowitz 
and W. Hill, published by Cambridge University Press, Cambridge, U.K. (1 989). Alternatively, the metrics may be meas- 
ured directly using analog electronics. The values of the metrics are calculated for a known segment of normal sinus 
rhythm (NSR). The mean of the NSR value (hereinafter "NSR-mean") for each metric is stored for later use. 

50 [0022] When a rhythm is to be diagnosed, the NSR-mean is subtracted from each measured metric in the unknown 
rhythm. This gives the NSR-mean-removed metrics. If the unknown rhythm is in fact ST, then the values of all the NSR- 
mean removed metrics will be small, because the shape of the electrogram will be similar to that in NSR. If the rhythm 
is VT, however, then at least one (and usually more) of the NSR-mean-removed metrics will be large. This difference 
is detected by the classification process, which is described below. 

55 [0023] The classification process includes a classifier, hereinafter described in greater detail and referred to as the 
"SUM classifier". The SUM classifier is an improvement over previously known methods of classification theory. It 
provides similar performance at a much reduced computational load compared to previously known methods. In a 
step-wise fashion, the SUM classifier performs the following operations: 
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a) Filters the electrogram; 

b) Calculates the raw metrics from the filtered electrogram; 

c) Amplitude-normalizes the metrics; 

d) Subtracts off the NSR-means; 
5 e) Subtracts off the ST-offsets; 

0 Takes the absolute values (to get error terms); 

g) Weights the error terms with weighing parameters; 

h) Sums the weighted error terms (to get the error sum); and 

i) Compares the error sum to a threshold. 

10 

The foregoing steps give a classification of successive heartbeats as either normal or abnormal. A normal rhythm (i. 
e., sinus tachycardia) is then defined as a rhythm that contains some number of normal heartbeats in a specified 
number of heartbeats, for example 7 or more normal heartbeats in the last 10 heartbeats. 

[0024] The metrics values are largely insensitive to filter corner frequencies. Sample rate also has little effect on 
15 performance, provided that peak-holding is employed. Good performance has been demonstrated at sample rates as 
low as 80 Hz, as indicated earlier. 

[0025] Quantization of the signal to 8 bits causes only a small degradation in performance, and this is the recom- 
mended signal resolution. It may be possible to quantize the signal further, but this has not been tested. 
[0026] The advantages of the device of the present invention for use in ICD's are, thus, that it employs a low sampling 
20 rate, requires very few multiplications per heartbeat, consumes low power, and can be implemented in analog elec- 
tronics. 

DESCRIPTION OF THE DRAWINGS 

25 [0027] While the specification concludes with claims particularly pointing out and distinctly claiming the subject matter 
of the invention, it is believed that the invention will be better understood from the following description, taken in con- 
junction with the accompanying drawings, in which: 

FIG. 1 shows the method of calculation of the MOD metric; 
so FIG. 2 shows the method of calculation of the PKP metric; 

FIG. 3 shows the method of calculation of the DIF metric; and 

FIGS. 4A-4I show the values of the QRS metric for a number of representative wave shapes. 
DESCRIPTION OF THE PREFERRED EMBODIMENT 

35 

[0028] A number of concepts recur in the following description of the metrics. Referring to FIG. 1, which illustrates 
a QRX complex 20 of an electrogram, the fiducial point or sensing point is shown at 10. The fiducial point 10 is a timing 
reference within the electrogram. In the simplest case the fiducial point is the point at which a sensing circuit (not 
shown) is tripped. The fiducial point that is currently used is the time of greatest dV/dt (slew) of the electrogram or, 

40 equivalent^, the time at which the high-pass-filtered electrogram has the greatest excursion from a baseline 12. 

[0029] An analysis window 14 is defined in terms of absolute offsets 16 and 18 from the fiducial point 10. In tests 
that have been performed utilizing the present method, the analysis window 14 has been set at 100 ms wide, centered 
on the fiducial point 10. The width of the analysis window 14 is not critical, but should extend out to the isoelectric 
segments on either side of the QRS complex 20 under consideration. Incoming (raw) data is filtered in an analysis filter 

45 (not shown) before any metric analysis is performed on the incoming data stream. The analysis filter removes low 
frequency base line wander, and high frequency noise. It also serves as an antialiasing filter. The analysis filter is 
preferably a first order 2 Hz highpass filter, followed by a second order 30 Hz Butterworth low pass filter. 
[0030] In the discussions which follow, references to an electrogram are intended to mean the electrogram at the 
output of the analysis filter, with the exception of the metric hereinafter referred to as TIL" metric, which metric uses 

so a bandpass filter (not shown, but hereinafter discussed) in place of the 30 Hz Butterworth low pass filter, 

[0031] The following metrics which may be used in connection with the present invention will now be described in 
greater detail: MOD, PK, PVE, NVE, PKP, DIF, FIL, SSA, QRS and ORD. 

MOD Metric 

55 

[0032] The MOD metric is the absolute (modulus) area under the QRS complex 20 of FIG. 1 . It is shown shaded at 
22 and 24 in that figure. The MOD metric is defined as the SUM of absolute values of the samples within the analysis 
window 14. For an analogue implementation, the MOD metric would be defined in terms of an integral. This metric is 
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intended to reflect the width of the complex 20 of the electrogram. 
PK Metric 

5 [0033] The PK metric is the value of the largest excursion of the electrogram complex 20 from the base line 1 2. It is 
signed so that if the largest peak is negative, than the PK metric is negative, and if the largest peak is positive (see, 
e.g. deflection 28 in FIG. 3), then the PK metric is positive. It reflects both the amplitude and polarity of the electrogram 
complex 20. 

10 PVE and NVE metrics 

[0034] The PVE metric is the amplitude of the largest positive deflection (see, e.g. deflection 28 in FIG. 3) in the 
electrogram complex 20. Similarly, the NVE metric is the amplitude of the largest negative deflection (see, e.g., de- 
flection 30 in FIG. 3) in the electrogram complex 20. 

15 

PKP metric 

[0035] The PKP metric is the peak-to-peak amplitude of the electrogram complex, and is shown in FIG. 2 at 26. This 
is the absolute amplitude of the largest positive deflection plus the absolute amplitude of the largest negative deflection. 
20 it reflects the amplitude of the electrogram complex 20. 

DIF metric 

[0036] Referring to FIG. 3, the DIF metric is the difference between the absolute amplitudes of the largest positive 
25 deflection 28 of the electrogram complex 20 and the largest negative deflection 30 of the electrogram complex 20. The 
DIF metric will be large and positive for an upright electrogram complex, small for a biphasic one, and large and negative 
for an inverted one. Thus, the DIF metric is a continuous variable which reflects the polarity of the QRS complex of the 
electrogram. 

30 FIL metric 

[0037] The FIL metric is the peak-to-peak amplitude of a filtered version of the electrogram. The raw signal is passed 
through a filter having the following characteristics: single pole high pass at 2 Hz, two pole band pass, and Q = 1 at 
25 Hz. The peak-to-peak amplitude of this signal is measured in the same way that the PKP metric is measured. This 
35 metric is designed to reflect the slew rate (dV/dt) of the electrogram complex 20. The high pass characteristic in the 
filter acts as a differentiator. 

SSA metric 

40 [0038] The SSA metric is the signed-square area metric. It is an integral measure, like the MOD metric. It is calculated 
as follows: for each sample within the analysis window the sample is squared and its sign is restored (e.g., 2 2 = 4; -2 2 
= -4). Then the signed squares are summed, the square root of the absolute value of the sum is taken, and the sign 
is restored. This gives a measure having characteristics similar to those of the DIF metric, i.e., it forms a continuous 
variable that reflects the polarity of the electrogram. The squaring is used to ensure that a large amplitude narrow lobe 

45 in the electrogram will predominate over a smaller amplitude but wider one. This is necessary because the high pass 
characteristic in the analysis filter will often produce a wide, low amplitude lobe in the filtered electrogram. The following 
is a mathematical description of the SSA metric: 

so t=S(sign(x).x 2 ) 

SSA = sign(t).7|t| 

55 QRS metric 

[0039] Referring to FIGS 4A-4I, the QRS metric represents the QRS angle / as calculated from the amplitude and 
temporal order of the largest positive and largest negative peaks, as follows: 
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r = abs(smaller peak + larger peak) 



H = arcos(r) 

[0040] if the first peak is positive and larger: /=H 
If the first peak is positive and smaller: /=-H 
If the first peak is negative and larger: /=H -180 
If the first peak is negative and smaller: /=180-H 

ORD metric 

[0041] The ORD metric reflects the order of the peaks in the electrogram complex 20. It is used in conjunction with 
the DIF metric to replace the QRS metric. The ORD metric is defined as follows: 
[0042] If the positive peak is first, ORD = smaller peak + larger peak. 
[0043] If the negative peak is first ORD = -smaller peak + larger peak. 

[0044] This gives a metric (ORD) which is a continuous variable that reflects order. Thus, if the electrogram complex 
20 is biphasic with the positive peak first, the result is close to one; if it is biphasic with the negative peak first, the result 
is close to minus one; if it is monophasic, the result is close to zero. 

[0045] The ORD and DIF metrics, in combination, would give ordered pairs of numbers that correspond very roughly 
to the QRS angles / of FIGS. 4A-4I, as shown in Table 1 below: 



Table 1 - 



Relationship Among QRS, ORD and DIF Metrics 


QRS angle / 


(DIF.ORD) 


0 


(0.0,1.0) 


45 


(0.5,0.5) 


90 


(1.0,0.0) 


135 


(0.5,-0.5) 


180 (=-180) 


(0.0,-1.0) 


-135 


(0.5,-0.5) 


-90 


(-1.0,0.0) 


-45 


(-0.5,0.5) 



[0046] All three of the metrics QRS, ORD and DIF suffer from the problem that they represent a two-point description 
of the electrogram complex 20, and some electrogram complexes necessarily require three points to properly describe 
them. They are simple to implement, however, requiring only positive and negative peak pickers (not shown), which 
are needed in any event for other metrics. 

[0047] As will be discussed in greater detail hereinafter, tests were performed to determine the correlation between 
various parameters. All of the metrics, except QRS angle /, were highly correlated with each other. This is because 
they all depend very strongly on the amplitude of the electrogram complex 20. This does not represent a problem for 
the PK and PKP metrics, which are intended to reflect the amplitude of the electrogram. However, it is unsatisfactory 
for the other metrics, which end up reflecting the amplitude more strongly than the desired feature. The result is a large 
collection of amplitude measures, and very little else. 

Amplitude Normalization 

[0048] To resolve the foregoing problem, metrics that have an unwanted sensitivity to amplitude are normalized in 
accordance with this invention. Such normalization is achieved by dividing the metrics in question by the PKP metric. 
This results in a set of metrics which are only loosely correlated with each other. Thus, they should reflect different 
aspects of the shape of the electrogram complex, as is desired. It is preferable to always use the MOD, PVE, NVE, 
DIF, FIL and SSA metrics in the normalized form. 
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[0049] This process of amplitude normalization is preferably integrated into the SUM classifier, described below, 
enabling the number of computer operations to be minimized. 

SUM classifier 

5 

[0050] The SUM classifier for the metrics involves a complex process which decides if a particular measurement 
resides in a region of the measurement space that is considered normal (to be classed as ST). If the measurement 
falls outside of this region, it will classed as VT. Using the metrics MOD, PKP, and NVE, the measurement space will 
be three dimensional, and the normal region will be a three dimensional space around the mean of the ST measure- 
10 ments. 

[0051 ] The goal of the SUM classifier is to minimize the computational cost of its operation. This is primarily dependent 
on the number of multiplications made per heartbeat. 

[0052] Given a three dimensional feature vector, for example the three measurements MOD, PKP and NVE, denoted 
m, p and n, respectively, and a set of weights (a, p, and y) the SUM classifier is defined as: 

15 



20 



ST 

a {m-m^l +p Ip-p^l +y \n-n ST \ ^ T 

VT 



(1) 



[0053] The above expression means that if the weighted sum of the distance of the measurement from the ST mean 
(ST offsets) is less than some threshold (T), the measurement will be classified as ST; otherwise, it would be classified 
25 as VT. Stated another way, it means that if the measurement is similar to the ST mean it will be classified as ST, and 
if it is dissimilar, it will be classified as VT. 

[0054] The foregoing discussion assumes that the measured metrics would be used directly. In fact, they must be 
processed somewhat before being used in the SUM classifier. Firstly, the m (MOD) and n (NVE) metrics need to be 
normalized by division by the p (PKP) metric, as follows: 

30 

--J (2) 

35 "' = j5 (3) 

[0055] The NSR-mean is then subtracted from the normalized metrics: 

A. m A, 

40 m" = m' -m NSR = - - m NSR (4) 

P" = P-Pnsr (5) 

45 A A , 

n "="'-n'NSR = y n 'NSR (6) 

A 

where m' 0TS , for example, is the mean of the normalized value of the MOD metric, calculated in NSR. 
[0056] The discriminant function is thus: 

ST 

55 L VT 

This must be reorganized firstly to remove divisions, and then to minimize the number of multiplies that.must be per- 
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10 



15 



20 



25 



30 



formed beat-by-beat: 

where a, g, y, m ST , p gT and n^. are parameters (and 



' ' . ST ' 

i l»-P I +p| |p- <i5*s*+p£) I *^a-P (iJ^flft) I • P < 8 ) 



T T T 

so need not be calculated beat-by-beat), and 



. .<») 



are calculated by background processes during NSR. This updating preferably occurs at a programmable interval in 
the range of 10 minutes to 1 hour. Thus, the number of multiplications required per complex during a rhythm classifi- 
cation is six. 

35 [0057] The terms a, p and y are calculated as the differences between the inverse of the square root of the diagonal 
terms of the class covariance matrices K, described in greater detail below. 

[0058] The classifier parameters for use in an implant_are calculated^from a set of metrics produced by the real 
implant hardware and software. The parameters a, p, y, m" p^ T and n^. are preferably programmable, but should 
not be available to the physician. This will allow the performance of the system to be tuned. The threshold parameter 
40 T (or an offset on the threshold) could be made available to the physician. This would allow a trade-off between sen- 
sitivity and specificity. 

[0059] In order to more clearly understand the term "covariance matrices K", referred to above, the following tutorial 
material with respect to random variables is provided. 

[0060] The expectation or mean of a random vector "x" is a vector "m" defined as: 

45 

m = E[x) = C„xp x (x)dx (12) 

where p y (x) is the probability density at x. This is taught by Therrien, C.W., in the book "Decision Estimation and 
50 Classification", published by John Wiley and Sons, New York, 1 989, (hereinafter referred to as Therrien"), at equation 
2.41, page 24. 

[0061] When training a classifier from a limited data set, this is estimated as: 
55 n "ri 
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where x jk is the j'th observation vector for class i; Nj is the number of observations of class i; and rrij is the estimate of 
the mean vector of class i. 1 
[0062] The covariance matrix "K" of a random vector V with expectation ¥ is defined as: 

K=E[(x-m)(x-m) J ] (14) 

where T is a matrix transposition operator (Therrien, page 49); and E is an expectation operator (Therrien, equa- 
tion 12, page 54). (Also, see Therrien, equation 4.42, page 55). 
[0063] When training a classifier from a limited data set, this is estimated as: 




(IS) 



where Kj is the estimate of the covariance matrix of class i. 
[0064] The SUM classifier relies on a specific discriminant function. Accordingly, the following discussion, relating 
to decision rules and discriminant functions, is presented at this time. 
[0065] The decision rule for for the two-class case takes the form: 

.0),' 

f(?) * T. (16) 
<*>2 

[0066] The interpretation of this expression is as follows. If some function f() of the observation vector y is less than 
some threshold T, then assign the class label <o A to it; otherwise assign the class label coj. 

[0067] The SUM classifier is derived from a Gaussian quadratic classifier, which is well known in the literature. Ac- 
cordingly, the following discussion relating to the Gaussian quadratic classifier is provided. 

[0068] The Gaussian quadratic classifier provides a Bayes optimal solution if the class conditional probability density 
functions are Gaussian. If the density functions are not Gaussian, then the probability of error is not minimised. The 
resulting classifier provides a decision surface which is best matched to the second moment statistics. (See, Therrien, 
pages 96,97). 

[0069] A general two-class quadratic classifier is defined as follows: 

h (y) ^ T Aj>+Jb T 5>+ c T ( 17 ) 

> 

(See Therrien, equation 6.6, page 96). 

[0070] This resulting classifier will be Gaussian quadratic if the classifier parameters are defined as follows: 

A = K?-Kl (18) 



b = 2(/r 2 1 m 2 -K' 1 1 m 1 ) 



(19) 
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< - m^K^rriy - K^ 1 m 2 + InL^ 

|K 2 



(20) 



5 [0071] The derivation of the SUM classifier, starting with the equation for the discriminant function of the Gaussian 
quadratic classifier, is set forth below. It develops the equations to calculate the parameters of the sum classifier, viz. 
alpha, beta, and gamma. 

[0072] The discriminant function for the SUM classifier can be derived from the Gaussian quadratic discriminant 
function, which can be written as: 

10 



h{y) = (y-m,) r K?(y-mi)- (y-m 2 )V 2 1 (y-m 2 ) + In 



(21) 



15 [0073] In practice this will be compared to a threshold (for a two-class problem): 



20 



h(y) K T 
> 



(22) 



[0074] Thus the log of the covariance determinants can be considered as part of the threshold value T, resulting in 
25 the simpler form: 



h'(y) = (y-m^Kliy-m,) - (y-m 2 ) T K2 1 (y.m 2 ) 
30 [0075] If we set the class 2 (VT) mean to be equal to the class 1 (ST) mean, we can make the substitution: 



(23) 



/ = y-m 1 (24) 

35 [0076] The discriminant function thus becomes: 

W = y^(<-K2V (25) 
40 [0077] If we set the off-diagonal terms for the covariance matrices to be zero, this becomes (for three features): 
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which reduces to: 
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(27) 



which reduces to: 



15 



(28) 



20 [0078] From the above we see that the discriminant function is a a weighted sum-of-square distance measure. (It 
actually measures the distances from the ST class mean, which intuitively makes sense.) 

[0079] In the special case where: < a' < a* and < a* this can be re-cast using a sum-of-absolute- 
values distance measure, giving: 



25 



30 



35 



40 



(29) 



where the coefficients a, p, and y of the SUM classifier are thus calculated from the diagonal terms of the class cov- 
ariance matrices as follows: 



P = 





J_ 








a a2 




_ J_ 








c b2 



(30) 



(31) 



(32) 



45 



In practice it has been found sufficient to approximate this as follows: 

« = JL.JL 

c a1 °a2 



(33) 
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y = JL.J_ 

G d <*c2 



(35) 



[0080] It will be apparent from the foregoing description that the present invention provides a device for classifying 
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heart rhythms by means of electrogram morphology, which device may be used in implantable cardioverter defibrillators, 
involves a low sample rate, very few multiplications per heart beat, can be implemented in digital as well as analog 
electronics, and consumes very little power. 

[0081] While particular embodiments of this invention have been shown and described, it will be obvious to those 
skilled in the art that various other changes and modifications may be made without departing from this invention. For 
example, the MOD, PKP and NVE metrics are not the only suitable set of metrics that may be employed with this 
invention. Many other combinations of the metrics described herein may also be used. Similarly, other metrics may be 
devised that might also work well. In addition, the parameters of the SUM classifier could be calculated using equations 
30-32, rather than 33-35. Alternatively, any other type of classifier known to the art could be substituted for the SUM 
classifier. These include discriminant function, K-nearest neighbors, neural network, ...ID3. It is, therefore, aimed in 
the appended claims to cover all such changes and modifications as fall within and scope of this invention. 



Claims 

1. A device for classifying heart rhythms of a patient by means of electrogram morphology, comprising: 

means for providing an electrogram of a patient's heart rhythm; 

means for determining at least two morphology-defining metrics derived from a QRS complex (20) of said 
electrogram, said electrogram having a width and an amplitude, said metrics defining the shape of said elec- 
trogram including at least one of said width and amplitude; and 

means for comparing the metrics of said given heartbeat with the corresponding metrics of a stored set of said 
at least two morphology-defining metrics of a normal electrogram; 

characterised by: 

means for establishing bounds of normal variation for said metrics which are combined to define a decision 
surface that includes first and second regions such that said first region represents normal heartbeats and 
said second region represents abnormal heartbeats, established from said stored set of said at least two 
morphology-defining metrics of a normal electrogram; 

wherein said comparing means is operable to decide which of said regions said given heartbeat falls into 
and to classify said given heartbeat as either normal or abnormal based on said decision; and 

means for diagnosing the patient's heart rhythm as ventricular arrhythmia if a number of said given heartbeats 
are classified as abnormal. 

2. A device according to claim 1 , wherein said determining means further includes means for determining a first set 
of at least three metrics selected from a group of metrics that includes: a MOD metric (22, 24), which comprises 
the absolute area under a QRS complex (20) of said electrogram; a PK metric (28), which comprises the value of 
the largest excursion of said QRS complex (20) from its baseline, said value having the same polarity as said 
excursion; a PVE metric (28), which comprises the amplitude of the largest positive deflection of said QRS complex 
(20); a NVE metric (30), which comprises the amplitude of the largest negative deflection of said QRS complex 
(20); a PKP metric, which comprises the peak-to-peak amplitude of said QRS complex (20); a DIP metric, which 
comprises the difference between the absolute amplitudes of the largest positive deflection of said QRS complex 
(20) and the largest negative deflection of said QRS complex (20); a FIL metric, which comprises the peak-to-peak 
amplitude of said QRS complex (20) after said complex has passed through a bandpass filter having a center 
frequency of about 25 Hz; a SSA metric, which comprises a signed-square area metric in which a.plurality of 
samples taken of said QRS complex (20) are squared, their signs are restored, the signed squares are summed, 
the square root of the absolute value of the sum is taken, and the sign is restored thereto; a QRS metric, which 
comprises a QRS angle © calculated from the amplitude and temporal order of the largest positive and largest 
negative peaks of the QRS complex (20), as follows: 

r = abs (smaller peak) + (larger peak), 



4> = arcos (r), 
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and 

if the first peak is positive and larger, 0 = 
if the first peak is positive and smaller, G = h|>, 
if the first peak is negative and larger, 6 = $ - 180, and, 
5 if the first peak is negative and smaller, 6 = 180 - <)>; 

and, an ORD metric, which is defined as follows: 

if the positive peak is first; ORD = smaller peak + larger peak, and 
if the negative peak is first, ORD = - smaller peak + larger peak. 

10 

3. A device according to claim 2, wherein said determining means includes means for determining a first set of at 
least three metrics that includes said MOD metric (22, 24), said PKP metric and said NVE metric (30) period. 

4. A device according to claim 2 t wherein, if one or more of the MOD, PVE, DIF, NVE, FIL and SSA metrics are 
15 selected from said group of metrics, said determining means is further operable for amplitude normalizing said 

one or more metrics by dividing the value thereof by the value of said PKP metric. 

5. A device according to any one of claims 1 to 4 further including means for filtering the patient's electrogram through 
a bandpass filter having a center frequency of about 25 Hz, and means for sampling said filtered electrogram at 

20 a predetermined rate to obtain said metric. 

6. A device according to claim 5 wherein said filtering means comprises a filter having a pass band of about 2 to 30 Hz. 

7. A device according to claim 6, wherein said sampling means includes means for storing a peak value of said 
25 electrogram, and wherein said predetermined rate is about 80 Hz. 

8. A device according to claim 6, wherein said determining means includes means for calculating and storing an 
NSR-mean value for each of the metrics of said normal electrogram, said NSR-mean value for each metric com- 
prising the mean value of the metric during normal sinus rhythm of a normal electrogram. 

30 

9. A device according to claim 8, wherein said comparing means includes means for subtracting from each of the 
metrics of said given heartbeat the corresponding NSR-mean value therefor to provide NSR-mean-removed met- 
rics corresponding to the patient's heart rhythm. 

35 10. A device according to claim 9, wherein said comparing means includes means for classifying a heartbeat as either 
normal or abnormal based on the value of one or more of said NSR-mean-removed metrics exceeding a prede- 
termined threshold value. 

11. A device according to claim 10, wherein said diagnosing means includes means for diagnosing the patient's heart 
<o rhythm as ventricular fibrillation if at least four or more heartbeats often successive heartbeats are classified as 

abnormal. 

12. A device according to any one of claims 1-11 , wherein said diagnosing means includes means for diagnosing the 
patient's heart rhythm as sinus tachycardia if at least a specified number of a predetermined greater number of 

45 last heartbeats are classified as normal. 

13. A device according to any one of claims 1-11, wherein said diagnosing means includes means for diagnosing the 
patient's heart rhythm as sinus tachycardia if at least 7 of the last 1 0 heartbeats are classified as normal. 



Patentanspruche 

1. Vorrichtung zum Klassifizieren von Herzrhythmen eines Patienten mittels einer Elektrogramm-Morphologie, mit: 

einer Einrichtung zur Bereitstellung eines Elektrogramms des Herzrhythmus eines Patienten; 
einer Einrichtung zur Bestimmung von mindestens zwei morphologiedefinierenden Metriken, die von einem 
QRS-Komplex (20) des Elektrogramms abgeleitet sind, wobei das Elektrogramm eine Breite und eine Ampli- 
tude aufweist, wobei die Metriken die Form des Elektrogramms definieren, das mindestens eines, namlich die 
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Breite und/oder die Amplitude aufweist; und 

einer Einrichtung zum Vergleichen der Metriken des gegebenen Herzschlags mit den entsprechenden Metri- 
ken eines gespeicherten Satzes der mindestens zwei morphologiedefinierenden Metriken eines normalen 
Elektrogramms; 

gekennzeichnet durch: 

eine Einrichtung zum Festlegen von Grenzen einer normalen Schwankung fur die Metriken, die kombiniert 
werden, urn eine Entscheidungsoberflache zu definieren, die einen ersten und einen zweiten Bereich aufweist, 
so daft der erste Bereich normale Herzschlage darstellt und der zweite Bereich unnormale Herzschlage dar- 
stellt, die anhand des gespeicherten Satzes der mindestens zwei morphologiedefinierenden Metriken eines 
normalen Elektrogramms festgelegt werden; 

wobei die Vergleichseinrichtung betriebsfahig ist zu entscheiden, in welchen von den Bereichen der gege- 
bene Herzschlag fallt, und den gegebenen Herzschlag auf der Grundlage der Entscheidung entweder als normal 
Oder als unnormal zu klassifizieren; und 

eine Einrichtung zum Diagnostizieren des Herzrhythmus eines Patienten als ventrikulare Arrhythmie, wenn 
eine Anzahl der gegebenen Herzschlage als unnormal klassifiziert wird. 

Vorrichtung nach Anspruch 1, wobei die Bestimmungseinrichtung ferner eine Einrichtung zum Bestimmen eines 
ersten Satzes von mindestens drei Metriken aufweist, die aus einer Gruppe von Metriken gewahlt sind, die aufweist: 
eine MOD-Metrik (22, 24), die die absolute Flache unter einem QRS-Komplex (20) des Elektrogramms umfaBt; 
eine PK-Metrik (28), die den Wert des groBten Ausschlags des QRS-Komplexes (20) von seiner Basislinie umfaBt, 
wobei der Wert die gleiche Polaritat hat wie der Ausschlag; eine PVE-Metrik (28), die die Amplitude mit der groBten 
positiven Ablenkung des QRS-Komplexes (20) umfaBt; eine NVE-Metrik (30), die die Amplitude der groBten ne- 
gativen Ablenkung des QRS-Komplexes (20) umfaBt; eine PKP-Metrik, die die Spitze-Spitze-Amplitude des 
QRS-Komplexes (20) umfaBt; eine DIF-Metrik, die die Differenz zwischen den absoluten Amplituden mit der 
groBten positiven Ablenkung des QRS-Komplexes (20) und der groBten negativen Ablenkung des QRS-Komple- 
xes (20) umfaBt; eine FIL-Metrik, die die Spitze-Spitze-Amplitude des QRS-Komplexes (20) umfaBt, nachdem der 
Komplex durch ein BandpaBfilter mit einer Mittenfrequenz von etwa 25 Hz gelaufen ist; eine SSA-Metrik, die eine 
vorzeichenbehaftete Quadratflachenmetrik umfaBt, bei der eine Vielzahl von Abtastwerten, die dem QRS-Komplex 
(20) entnommen werden, quadriert werden, deren Vorzeichen wiederhergestellt werden, die vorzeichenbehafteten 
Quadrate summiert werden, die Quadratwurzel des absoluten Wertes der Summe gebildet wird und das Vorzeichen 
dazu wiederhergestellt wird; eine QRS-Metrik, die einen QRS-Winkel e umfaBt, der aus der Amplitude und der 
zeitlichen Reihenfolge der groBten positiven und der groBten negativen Spitze des QRS-Komplexes (20) berechnet 
wird, namlich wie folgt: 

r = abs (kleinere Spitze) + (groBere Spitze), 



<t> = arc (r), 

und 

wenn die erste Spitze positiv und groBer ist, 0 = 

wenn die erste Spitze positiv und kleiner ist, 0 = 

wenn die erste Spitze negativ und groBer ist, 0 = <(> - 180 und 

wenn die erste Spitze negativ und kleiner ist, 0 = 180 - (J>; 

und eine ORD-Metrik, die folgendermaBen definiert ist: 

Wenn die positive Spitze zuerst auftritt, ORD = (kleinere Spitze) + (groBere Spitze), und 
wenn die negative Spitze zuerst auftritt, ORD = -(kleinere Spitze) + (groBere Spitze). 

Vorrichtung nach Anspruch 2, wobei die Bestimmungseinrichtung eine Einrichtung zum Bestimmen eines ersten 
Satzes von mindestens drei Metriken aufweist, die die Periode der MOD-Metrik (22, 24), der PKP-Metrik und der 
NVE-Metrik (30) aufweist. 

Vorrichtung nach Anspruch 2, wobei, wenn eine oder mehrere der MOD-, PVE-, DIF-, NVE-, FIL- und SSA-Metriken 
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aus der Gruppe von Metriken gewahlt werden, die Bestimmungseinrichtung ferner fur eine Amplitudennormierung 
der einen oder mehreren Metriken durch Teilung ihres Wertes durch den Wert der PKP-Metrik betriebsfahig ist. 

5. Vorrichtung nach einem der Anspruche 1 bis 4, ferner mit einer Einrichtung zum Filtern eines Elektrogramms des 
5 Patienten mittels eines BandpaBfilters mit einer Mittenfrequenz von etwa 25 Hz und einer Einrichtung zum Abtasten 

des gefilterten Elektrogramms mit einer vorbestimmten Rate, um die Metrik zu ermitteln. 

6. Vorrichtung nach Anspruch 5, wobei die Filtereinrichtung ein Filter mit einem BandpaB von etwa 2 bis 30 Hz umfaBt. 

10 7. Vorrichtung nach Anspruch 6, wobei die Abtasteinrichtung eine Einrichtung zum Speichern eines Spitzenwerts 
des Elektrogramms aufweist und wobei die vorbestimmte Rate etwa 80 Hz ist. 

8. Vorrichtung nach Anspruch 6, wobei die Bestimmungseinrichtung eine Einrichtung zum Berechnen und Speichern 
eines NSR-Mittelwertes fur jede der Metriken des normalen Elektrogramms aufweist, wobei der NSR-Mittelwert 

is fur jede Metrik den Mittelwert der Metrik wahrend eines normalen Sinusrhythmus eines normalen Elektrogramms 

umfa&t. 

9. Vorrichtung nach Anspruch 8, wobei die Vergleichseinrichtung aufweist: eine Einrichtung zum Subtrahieren des 
entsprechenden NSR-Mittelwertes des gegebenen Herzschlags von jeder der Metriken, um NSR-Mittelwert-be- 

20 reinigte Metriken bereitzustellen, die dem Herzrhythmus des Patienten entsprechen. 

10. Vorrichtung nach Anspruch 9, wobei die Vergleichseinrichtung aufweist: eine Einrichtung zum Klassifizieren eines 
Herzschlags entweder als normal oder unnormal auf der Grundlage des Wertes einer oder mehrerer der NSR-Mit- 
telwert-bereinigten Metriken, der einen vorbestimmten Schwellwert uberschreitet. 

25 

1 1 . Vorrichtung nach Anspruch 1 0, wobei die Diagnoseeinrichtung aufweist: eine Einrichtung zum Diagnostizieren des 
Herzrhythmus eines Patienten als Harzkammerflimmern, wenn mindestens vieroder mehr Herzschlage von zehn 
aufeinanderfolgenden Herzschlagen als unnormal klassifiziert werden. 

30 12. Vorrichtung nach einem der Anspruche 1 bis 11, wobei die Diagnoseeinrichtung aufweist: eine Einrichtung zum 
Diagnostizieren des Herzrhythmus eines Patienten als Sinustachykardie, wenn mindestens eine vorgegebene 
Anzahl einer vorbestimmten groBeren Anzahl der letzten Herzschlage als normal klassifiziert wird. 

13. Vorrichtung nach einem der Anspruche 1 bis 11, wobei die Diagnoseeinrichtung aufweist: eine Einrichtung zum 
35 Diagnostizieren des Herzrhythmus eines Patienten als Sinustachykardie, wenn mindestens sieben der letzten 

zehn Herzschlage als normal klassifiziert werden. 



Revendications 

40 

1 . Dispositif pour classer des rythmes cardiaques d'un patient au moyen d'une morphologie d'electrocardiogramme, 
comprenant: 
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un moyen pour produire un Electrocardiogramme d'un rythme cardiaque de patient; 

un moyen pour determiner au moins deux mesures de definition de morphologie qui sont deriv6es d partir 
d'un complexe QRS (20) dudit Electrocardiogramme, ledit electrocardiogramme presentant une largeur et une 
amplitude, lesdites mesures definissant la forme dudit electrocardiogramme incluant au moins un element 
pris parmi ladite largeur et ladite amplitude; et 

un moyen pour comparer les mesures dudit battement cardiaque donne aux mesures correspondantes d'un 
jeu stocke desdites au moins deux mesures de definition de morphologie d'un electrocardiogramme normal, 

caracterise par 

un moyen pour etablir des fimites d'une variation normale pour lesdites mesures, qui sont combiners pour 
definir une surface de decision qui inclut des premiere et seconde regions de telle sorte que ladite premiere 
region represente des battements cardiaques normaux et que ladite seconde region represente des batte- 
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ments cardiaques anormaux, comme etabli a partir dudit jeu stocke desdites au moins deux mesures de 
definition de morphologie d'un electrocardiogramme normal, 

dans lequel tedit moyen de comparaison peut fonctionner pour decider a i'interieur de laquelle desdites 
regions ledit battement cardiaque donne tombe et pour classer ledit battement cardiaque donne soit en tant que 
battement cardiaque normal, soit en tant que battement cardiaque anormal sur la base de ladite decision; et 

un moyen pour diagnostiquer le rythme cardiaque du patient en tant qu'arythmie ventriculaire si un certain 
nombre desdits battements cardiaques donnes sont classes comme etant anormaux. 

Dispositif selon la revendication 1, dans lequel ledit moyen de determination inclut en outre un moyen pour deter- 
miner un premier jeu d'au moins trois mesures choisies parmi un groupe de mesures qui inclut: une mesure MOD 
(22, 24) qui comprend I'aire absolue sous un complexe QRS (20) dudit electrocardiogramme; une mesure PK (28) 
qui comprend la valeur de I'excursion la plus grande dudit complexe QRS (20) par rapport a sa ligne de base, 
ladite valeur pr£sentant la m§me polarite que ladite excursion; une mesure PVE (28) qui comprend Pamplitude de 
la deviation positive la plus grande dudit complexe QRS (20); une mesure NVE (30) qui comprend ramplitude de 
la deviation negative la plus grande dudit complexe QRS (20); une mesure PKP qui comprend Pamplitude crete 
a cr§te dudit complexe QRS (20); une mesure DIF qui comprend la difference entre les amplitudes absolues de 
la deviation positive la plus grande dudit complexe QRS (20) et de la deviation negative la plus grande dudit 
complexe QRS (20); une mesure FIL qui comprend Pamplitude cr§te a crete dudit complexe QRS (20) aprfcs que 
ledit complexe est pass6 au travers d'un filtre passe-bande qui presente une frequence centrale d'environ 25 Hz; 
une mesure SSA qui comprend une mesure d'aire au carre sign6e ou une plurality d'echantillons extraits a partir 
du complexe QRS (20) sont eieves au carr6, leurs signes sont restaurs, les carr6s signes sont sommes, la racine 
carre de la valeur absolue de la somme est prise et le signe est restaure dessus; une mesure QRS qui comprend 
un angle QRS 0 qui est calcuie a partir de Pamplitude et de I'ordre temporel des cretes positives la plus grande 
et negative la plus grande du complexe QRS (20) comme suit: 

r = abs (crete plus petite)/(crete plus grande) 



$ = arcos (r); 

et 

si la premiere crete est positive et plus grande, e = O, 

si la premiere crete est positive et plus petite 6 = - <D, 

si la premiere crete est negative et plus grande, 8 = O - 180; et 

si la premiere crete est negative et plus petite, 6 = 180 - O; et 

une mesure ORD qui est d6finie comme suit: 

si la crete positive est en premier, ORD = crete plus petite/crete plus grande; et 

si la crete negative est en premier, ORD = - crete plus petite/crete plus grande. 

Dispositif selon la revendication 2, dans lequel ledit moyen de determination inclut un moyen pour determiner un 
premier jeu d'au moins trois mesures qui incluent ladite mesure MOD (22, 24), ladite mesure PKP et ladite periode 
de mesure NVE (30). 

Dispositif selon la revendication 2, dans lequel, si une ou plusieurs des mesures MOD, PVE, DIF, NVE, FIL et SSA 
sont choisies parmi ledit groupe de mesures, ledit moyen de determination peut en outre fonctionner pour norma- 
liser en amplitude lesdites une ou plusieurs mesures en divisant sa valeur ou leurs valeurs par la valeur de ladite 
mesure PKP. 

Dispositif selon I'une quelconque des revendications 1 a 4, incluant en outre un moyen pour filtrer Peiectrocardio- 
gramme de patient au travers d'un filtre passe-bande qui presente une frequence centrale a environ 25 Hz, et un 
moyen pour 6chantillonner ledit electrocardiogramme filtre a une frequence predetermin6e pour obtenir ladite 
mesure. 

Dispositif selon la revendication 5, dans lequel ledit moyen de filtrage comprend un filtre qui presente une bande 
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passante d'environ 2 a 30 Hz. 

7. Dispositif selon la revendication 6, dans lequel ledit moyen d'echantillonnage inclut un moyen pour stacker une 
valeur de crete dudit electrocardiogramme et dans lequel ladite frequence predeterminee est d'environ 80 Hz. 

8. Dispositif selon la revendication 6, dans lequel ledit moyen de determination inclut un moyen pour calculer et 
stocker une valeur moyenne NSR pour chacune des mesures dudit electrocardiogramme normal, ladite valeur 
moyenne NSR pour chaque mesure comprenant la valeur moyenne de la mesure pendant un rythme sinusal 
normal d'un electrocardiogramme normal. 

9. Dispositif selon la revendication 8, dans lequel ledit moyen de comparaison inclut un moyen pour soustraire de 
chacune des mesures dudit battement cardiaque donne la valeur moyenne NSR correspondante afferente afin de 
produire des mesures a moyenne NSR en levee correspondant au rythme cardiaque du patient. 

15 10. Dispositif selon la revendication 9, dans lequel ledit moyen de comparaison inclut un moyen pour classer un bat- 
tement cardiaque soit en tant que battement cardiaque normal, soit en tant que battement cardiaque anormal sur 
la base de la valeur d'une ou de plusieurs desdites valeurs a moyenne NSR enlevee excedant une valeur de seuil 
predeterminee. 

20 11. Dispositif selon la revendication 10, dans lequel ledit moyen de diagnostic inclut un moyen pour diagnostiquer le 
rythme cardiaque du patient en tant que fibrillation ventriculaire si au moins quatre battements cardiaques ou plus 
pris parmi dix battements cardiaques successifs sont classes comme etant anormaux. 

12. Dispositif selon Tune quelconque des revendications 1 a 11 , dans lequel ledit moyen de diagnostic Inclut un moyen 
25 pour diagnostiquer le rythme cardiaque du patient en tant que tachycardie sinusale si au moins un nombre specifie 

d'un nombre plus grand predetermine de derniers battements cardiaques sont classes comme etant normaux. 

13. Dispositif selon Tune quelconque des revendications 1 a 11, dans lequel ledit moyen de diagnostic inclut un moyen 
pour diagnostiquer le rythme cardiaque du patient en tant que tachycardie sinusale si au moins sept des dix derniers 

30 battements cardiaques sont classes en tant que battements cardiaques normaux. 
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